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Abstract 

Nonlinear contraction theory is a comparatively recent dynamic 
control system design tool based on an exact differential analysis of 
convergence, in essence converting a nonlinear stability problem into 
a linear time-varying stability problem. Contraction analysis relies 
on finding a suitable metric to study a generally nonlinear and time- 
varying system. This paper shows that the computation of the metric 
may be largely simplified or indeed avoided altogether by extending 
the exact differential analysis to the higher-order dynamics of the non- 
linear system. Simple applications in economics, classical mechanics, 
and process control are described. 



1 Introduction 



Nonlinear contraction theory is a comparatively recent dynamic control sys- 
tem design tool based on an exact differential analysis of convergence [T2"] . 
In essence, it allows one to convert a nonlinear stability problem into a lin- 
ear time-varying stability problem. Contraction analysis relies on finding a 
suitable metric to study a generally nonlinear and time-varying system. De- 
pending on the application, the metric may be trivial (identity or rescaling of 
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states), or obtained from physics, combination of contracting subsystems [T2] . 
semi-definite programming fX3j, or recently sums-of-squares programming [2J. 

The goal of this paper is to show that the computation of the metric may 
be largely simplified or avoided altogether by extending the exact differential 
analysis to the higher-order dynamics of the nonlinear system. Intuitively 
this is not surprising, since, as an elementary instance, a scalar linear time- 
invariant system would require in the original approach a non-identity metric 
(obtained from a Lyapunov Matrix Equation). 

After a brief review of contraction theory in Section 2, the main results 
are presented in Section 3, first in the discrete-time case (with a simple 
application to price dynamics in economics) and then in the continuous-time 
case. Simple examples and applications are discussed in Section 4, in the 
contexts of classical mechanics, process control, and observer design (see 
PJ El EE 120] for other recent applications of contraction theory to observer 
design). Hamiltonian systems are studied in section 5. Concluding remarks 
are offered in section 6. 



2 Contraction theory 

The basic theorem of contraction analysis ^2] can be stated as 

Theorem 1 Consider the deterministic system x = f(x, t) , where f is a 
smooth nonlinear function. If there exist a uniformly positive definite metric 



such that the Hermitian part of the associated generalized Jacobian F = 



is uniformly negative definite, then all system trajectories converge exponen- 
tially to a single trajectory, with convergence rate \X max \, where X max is the 
largest eigenvalue of the Hermitian part of F. The system is said to be 
contracting. 

In the above, * denotes complex conjugation, ~* for a matrix denotes 
the inverse of the conjugate matrix, and the state-space is R n (in this paper) 



M(x,t) 



0(x,t) 



r e*(x,t) 
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or C n . The system is said to be semi- contracting (for the metric M(x, t)) if 
F is always negative semi-definite, and indifferent if F is always zero. 

It can be shown conversely that the existence of a uniformly positive 
definite metric with respect to which the system is contracting is also a 
necessary condition for global exponential convergence of trajectories. In the 
linear time-invariant system is globally contracting if and only if it 

is strictly stable, with F simply being a normal Jordan form of the system 
and © the coordinate transformation to that form. Conceptually, approaches 
closely related to contraction, although not based on differential analysis, can 
be traced back to jS] and even to [TTj . 

Similarly, a discrete system 

x i+l = fi( x ij i) 

will be contracting in a metric T @* if the largest singular value of the 
discrete Jacobian 0j + i^-@~* is strictly smaller than 1. In the particular 
case of real autonomous systems with identity metric, the basic contraction 
theorem corresponds in the continuous-time case to Krasovkii's theorem [T5] . 
and in the discrete-time case to the contraction mapping theorem 

Contraction theory proofs make extensive use of virtual displacements, 
which are differential displacements at fixed time borrowed from mathemat- 
ical physics and optimization theory. Formally, if we view the position of 
the system at time t as a smooth function of the initial condition x and of 
time, x = x(x Q ,£) , then 5x = jj^- dx . For instance ^2] ; for the system 
of Theorem 1, one easily computes 

±(etx) = F(05x) (1) 

An appropriate metric to show that the system is contracting may be 
obtained from physics, combination of contracting subsystems ^2], semi- 
definite programming or sums-of-squares programming The goal of 
this paper is to show that the computation of the metric may be largely sim- 
plified or avoided altogether by considering the system's higher- order virtual 
dynamics (rather than merely its first-order virtual dynamics, as in equation 
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3 Higher-order contraction 

3.1 The discrete-time case 

Technically, the extension to higher-order contraction is simplest in the discrete- 
time case, which we discuss first. The main idea is to construct an exponential 
bound on the virtual displacement 5x over n time-steps, rather than over a 
single time-step as in |12j . 

Consider for i > the 72,-dimensional (n > 1) virtual dynamics 
5x i+n = A" -1 5x i+n _ 1 + . . . + A°5 Xi 
Taking the norm (denoted by | | ) on both sides, and bounding, yields 

\S*i+n\ < l^r 1 ! \ fei+n-l\ + ••• + l<N 

where the norm of a matrix is the largest singular value of that matrix. Let 
us bound for i = the initial conditions using real positive constants A and 
K as 

\6xj\ < KX j , 0<j<n 
Assume now that the following characteristic equation is verified, 

A n > |A"- 1 |A n - 1 + . . . + V? > 

We then get 

|5x i+n | < lA^lKX 1 ^- 1 + ... + \A° i \KX i < K\ l+n 
Repeating the above recursively for i > we get by complete induction 

|<5xi| < KX 

and hence exponential convergence of |<5xj|, as illustrated in Figure El 

Theorem 2 Consider j or % > the n- dimensional (n > 1) virtual dynamics 
5x i+n = A™- 1 ^^.! + . . . + A?5xi 
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6xi 








Figure 1: ±XA* defined by 5xi over i 

Let us define Vi > a constant X with the characteristic equation 

X n > \A^ 1 \X n - 1 + ... + \A°\, Vi > (2) 
VFe can i/ien conclude 

|<5x i+n | < KX l+n 

where K is defined by 

\5*j\ < KX j , < j < n (3) 
Thus, the system is contracting if X < 1. 

Example 3.1: Consider first a second-order linear time invariant (LTI) dy- 
namics 

X l+2 + 27 X i+1 + 07 2 Xj = Ui 

where Ui is an input, and 7 and o are constants. The virtual dynamics is 

5x i+2 = -27 Sxi+i - «7 2 5xi 
The characteristic equation (j2j) for A > is then given by 

A 2 > 2|7|A + |a|7 2 i.e. A > |7|(1 + v 7 ! + H ) 
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Thus, the contraction condition A < 1, or 



| 7 |(1 + y/l + \a\ ) < 1 

simply means that both eigenvalues of the system have to lie for the conjugate 
complex case (a > 1) within the red half circles in @ or on the green line for the 
real case (a < 1). Note that Theorem O simply bounds the possibly oscillating 
discrete system with a non-oscillating system of the same convergence rate for 
the real case. 

Consider now the virtual dynamics of an arbitrary second-order nonlinear time- 
varying system, 

5x i+2 + 2 7(i) 5x i+ i + a(i) j 2 (i)5xi = 

The characteristic equation and the contraction condition are the same as above, 
except that 7 and a are now time-dependent. D 




Figure 2: Contraction region in the complex plane of second-order discrete 
system 
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Example 3.2: In economics, consider the price dynamics 

n i+1 = fi{Pi,i) 
Pi+i = Ei( n i,i) 

with rij the number of sold products at time % and corresponding price pj . 

The first line above defines the customer demand as a reaction to a given price. 
The second line defines the price, given by the production cost under competi- 
tion, as a reaction to the number of sold items. The dynamics above corresponds 
to the second-order economic growth cycle dynamics 

n i+2 = fj+i (gi(xii,i)) 

Contraction behavior of this economic behavior with contraction rate A can then 
be concluded with Theorem |21 for 

That means we get stable (contraction) behavior if the product of customer 
demand sensitivity to price and production cost sensitivity to number of sold 
items has singular values less than 1. We can get unstable (diverging) behavior 
for the opposite case. 

Note that this result even holds when no precise model of the sensitivity is 
known, which is usually the case in economic or game situations. 

Whereas the above is well known for LTI economic models we can see that 
the economic behavior is unchanged for a non-linear, time-varying economic 
environment. 

Let as assume now that the above corresponds to a game situation (see e.g. 
[T%] or [3]) between two players with strategic action p^ and rij. Both players 
optimize their reaction g and f with respect to the opponent's action. 

We can then again conclude for Q to global contraction behavior to a unique 
time-dependent trajectory (in the autonomous case, the Nash equilibrium). D 

Of course, and throughout this paper, in some cases the analysis may yet 
be further streamlined by first applying a simplifying metric transformation 
of the form 5z = ©5x , and then applying the results to 5z. 
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3.2 The continuous-time case 



Let us now derive the continuous-time version of the previous results. Con- 
sider for t > the n-dimensional (n > 1) virtual dynamics 

tfx< ft) = -A n _ 1 5x^- 1 ) - ... - A 5x 

The proof is based on splitting up the dynamics into a stable part, described 
by a block diagonal matrix composed of identical negative definite blocks F 
which we select, and an unstable higher-order part. Let <5x Q = <5x, and define 
recursively 

<5x = F5x G + 5x1 



5± n ^ 2 = F5x n _ 2 + <5x„_i 

5± n ^ = -A n _Mr l) - ^n-2^ { r 2) - ••• -A 5x 
-(F5x ) ( - 1) -...-(F5x„_ 2 ) (1) 
= -A„_ 1 (F5x + <5x 1 ) (n - 2) -...-A 5x 

- (L l F5x + L°F5 Xl ) (n ~ 2) - ... - (L x F5x n _ 2 + m^) 
= F5x„_i - A;_ l( 5x n _i - A^_ 2 5x n _ 2 - A;_ 3 <5x n _3 - ... (5) 



where 



L°F 



L J F = — U~ l F + L j - l F F / > I 

dt 



and 



= A n _ 1+ ^jL°F (6) 
A;_ 2 = A n _ 2 + ( n ~ 1 ) A n _iL°F +f^L 1 F 

K_ 3 = A n _ 3 +( n 7 2 )A«- 2 L°F+( ?2 - 1 )A- 1 L 1 F+(ML 2 F 



Equation (jSJ) represents the superposition of a higher-order-system and a 
block diagonal dynamics in the chosen F. Let us assess the contraction 
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behavior of the higher-order part by taking the norm 

|*xW| < ia^ii^- 1 )] + |a;„ 2 ||^- 2 )| + . . . 

where the norm of a matrix is the largest singular value of that matrix. Let 
us bound for t = the initial conditions with real and constant A, K > and 
assume the following characteristic equation 

15x^1 < KX j e xt , 0<j<n (7) 

\ n > \K-i\* n ~ 1 + ■■■ + \Kl Vt>0 (8) 

Figure El shows how K has to be selected for a given A for a second-order 
system (n = 2). With (JHJ) we can bound the ro'th derivative of 5x as 

|5x (n) | < \A* n _ 1 \K\ n - l e xt + . . . + |A*|i^e A * < KX n e xt 

Integrating the above for t > we can exponentially bound the higher-order 
dynamics 5x as 

\8x.\ < Ke J ° XdT 
Using the above this allow to conclude: 

Theorem 3 Consider for t > the n- dimensional (n>l) virtual dynamics 

5x (n) = -A"- 1 ^"-" - ... - A°5x 

Let us define a constant A > such that Vt > we fulfill the characteristic 
equation 

A" > A^A- 1 + . . . + A: (9) 
where A* is defined in (OJ) for a given choice of the matrix F. 

We can then conclude on contraction rate (i.e., the largest eigenvalue of 
the symmetric part of) F + AI ; where |5x| is initially bounded with K , defined 
in 0. 

One specific choice of F is — A ^ -1 , which cancels the highest time-derivative 
on the right-hand side, and is known for LTV systems as the reduced or 
unstable form JU] of the original higher-order dynamics. We will use this 
definition of F in most of the following examples. Also note that more general 
forms could be chosen for the stable part. 
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Sx(t) 
. t 



Figure 3: Ke bounding the first and second derivative of 5x for a given A 
at t = 

4 Examples and Applications 

In this section, we discuss simple examples (section 4.1), applications to 
nonlinear observer design (section 4.2), and adding an indifferent system 
(section 4.3). 

4.1 Some simple examples 

Example 4.1: Consider the second-order LTI dynamics 

x = —2Qujx — oj 2 x 

with constant Q and uj > 0. The virtual dynamics is 

5x = —2(u)6x — u) 2 6x 

The characteristic equation Q is then given with F = Quj for constant, positive 
A by 
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Using Theorem |3] we can then conclude on contraction behavior with conver- 
gence rate 

This means that we require the poles to lie within the ±45° quadrant of the 
left-half complex plane. D 

While Theorem El can thus be overly conservative for LTI systems, this is 
not the case for general nonlinear time- varying systems, as we now illustrate. 

Example 4.2: Consider the second-order LTV dynamics 

x + a\ x + a Q {t) x = u(t) 

with ai,a > 0, which would be sufficient conditions for LTI stability. Let us 
assume a small damping gain a\ and strong spring gains a such that the system 
oscillates. 

If now the time- varying gain a a (t) is chosen to be very large when the system 
oscillates back to and small otherwise then the energy is constantly increased, 
which makes the system unstable (Figure EJ. This is precisely what is excluded 
by Theorem |31 d 




Figure 4: LTV system with a = 900 + 850 sin t and a\ = 0.01 
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Example 4.3: Consider the second-order LTV dynamics 

x + a\(t) x + a (t) x = u(t) 

The virtual dynamics is 

5x = —a\5x — a a 5x 
The characteristic equation ® is then given with F = — ^ by 

2 

. o i (Xi Cb\ , 

A 2 > + ^ + y| 

for constant positive A. Using Theorem El we can then conclude on contraction 
behavior with convergence rate + A (Figure EJ). □ 



POSITION AND VELOCITY 




Figure 5: LTV system with ao = 2 + cost and a\ = 8 + sint 



Example 4.4: Let us illustrate a case where choosing F other than — A ~ n ~ 1 
can simplify the result. Consider the generalized Van der Pol or Lienard dy- 
namics 

x = —a\{x)x — a a (x,t) 
with ati, > 0. The virtual dynamics is 

5x = —a\5x — ( -—^ + a'i ) 5x 
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The characteristic equation @ is then given with F = —a% + by 

2 . min(ai) 2 aimin(oi) <9a 
A = (ai - mm(a 1 ))A + | + — | 



and hence A > 0.5 ya\ — min(a 1 ) + y (ai — min(ai)) 2 + | — min(ai) 2 + 2a\ min(ai) — 4^| ) . 
Thus the contraction behavior of 5x is then given by 



da 

A + F < -0.5ai + 0.5W (ai - min(ai)) 2 + | - (01 - min(ai)) 2 + af - 4-rp| 

which is negative if the argument of the square root is less than af, i.e. for 
min(oi) 2 - 2ai min(ai) + 4^ < 0. 

Hence we can conclude on contraction behavior if the poles of the minimal 
damping min(ai) with the actual spring gain lie within the ±45° quadrant 
of the left-half complex plane (Figure EJ). Note that this result can be extended 
to the vector case x if the corresponding matrix -Ai(x) is integrable. D 



POSITIONS AND VELOCITIES 




-4 1 1 1 1 1 1 1 1 1 1 1 

O 10 20 30 40 50 60 70 80 90 100 

TIME(SEC) 



Figure 6: Van der Pol with a\ = 8 + sinx and u(t) = cost 

Example 4.5: Consider the second-order nonlinear vector system 

dV 

x + D x + — = u(t) 
ax 
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with potential energy V = x\ + x\ + Xi:r 2 sini and constant damping gain 
D = diag(l,4). 

The corresponding variational dynamics is 

d 2 v 

5x + D 5x + — — 7t 5x = 
ox 2 

The characteristic equation Q for F = — 5 is then given by 



2 

A > 



2 . ,3V DD 



3x 2 4 1 

for constant positive A. Using Theorem |3] we can then conclude on contraction 
behavior with convergence rate AI — S in Figure [7| D 



POSITIONS AND VELOCITIES 



01 23456789 10 

TIME(SEC) 

CONVERGENCE RATE 

-0.5 1 1 1 1 1 1 1 1 1 1 1 




01 23456789 10 

TIME(SEC) 



Figure 7: Second-order system with D = diag(l, 4) and V = x\ + x\ + x\X2 sin£ 

Example 4.6: Let us consider two coupled systems of the same dimensions, 
with the virtual dynamics 

d_ ( ci Xl \ / F n F 12 

dt V Sx 2 J V F 2i F 22 J V 5x 2 J 

Let us transform this dynamics in the following second-order dynamics 

d 2 d 

—2 <5xi = Fn<5xi + Fn — 5xi + (F12 + Fi 2 F 22 )5x 2 + Fi 2 F 2 i<5xi 

dt z dt 

= (Fn + F* 22 ) j t 5x 1 + (Fu + Fi 2 F 2 i - F* 22 F n ) <5xi 
14 



with the generalized Jacobian F22 = (F12 + F 12^22)1^ . The shrinking rate 
of this system is now the average of Fn and F| 2 . Using a direct contraction 
approach with e.g. F12 = — k F^ the guaranteed contraction rate would be a 
more conservative value, namely the largest of the individual contraction rates 
of Fn and F^. □ 



4.2 Higher-order observer design 

While a controller for an n th order system simply has to add stabilizing feed- 
back in x( n_1 \ ... , x, t according to Theorem El the situation is not such 
straightforward for observers since here only a part of the state is measured. 
Motivated by the linear Luenberger observer and the linear reduced-order 
Luenberger observer, we derive such an observer design for higher-order non- 
linear systems. 

Consider the n-dimensional nonlinear system dynamics 

xW = a (x,x,t)+ ... +air 2 2) (x,x,t) + ai'r i 1) (x,t) 

with the measurement y(x( n_1 ), ... , x, t). Note that for a linear Luenberger 
observer y is equivalent to x and all a^ are linear functions of x. 

Consider now the corresponding nonlinear observer 

x„-i = a - e (y) + e (y) 

x n -2 = x n _i + ai - ei(y) +ei(y) 

x = xi + a„_i - e n _i(y) + e n _i(y) 

with aj(x ,X! + a n _i — e„_i(y) + e„_i(y),t) and x G = x. Note that the 
coordinate transformation in the bracket is a nonlinear generalization of the 
reduced Luenberger observer. The above dynamics is equivalent to 

x<"> = a (x,x,t)+ ... +at- 2 2) (x,x,t) + a n n r i 1) (x,t) 
- e (y) + e (y)- ... - (y) + (y) 

whose variational dynamics is 

<5x^ = 5a + ... + Sa.^ -8e - ... - Se^ 
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where the varation is performed on x" , ... ,x. Hence the feedback is not 
only performed in y, but implicitly also up to the (n — l) th time-derivative 
of y. 



Example 4.7: Consider the second-order nonlinear system 

x + — — (x) x + a (x) = 
ox 

where x is measured. Consider now the corresponding nonlinear observer 

&i = a (x) — e Q (x — x) 

x = x\ + a\(x) — e\{x — x) 

with constant e Q and e\ and where we have replaced with the feedback a\ as a 
function of x. The corresponding second-order variational dynamics is 

da 

5x + ei 5x + (e Q — —r-) Sx = 
ox 

Contraction behavior can then be shown with Theorem |3 d 



Example 4.8: Consider the temperature-dependent reaction A — > B in a 
closed tank 

d / C,4 \ / — 1 \ _E 



dt\T J V -10 ' 6 TCA 

with ca the concentration of A, T the measured temperature, and E the specific 
activation energy. This reaction dynamics is equivalent to the following second- 
order dynamics in temperature 

— E ■ o E ■ 

T +YT T 2 = -e-rT 

rp — E 

Letting r = J o e~dT yields 

f = — e Tj- 

Contraction can then be shown as in Example 14.41 The observer dynamics 



i 



_ E_ X E ~n f T -E 

-e t oTo + —T% -e eTo / e ~ dT 
T„ 2 Jt 



T = Ti-ei(T -T ) 
with T Q = T and constant e Q and e\ leads to 

_ E 

f = (— e t - ei)f — e a f + e\f(t) + e r(i) 
whose contraction behavior can now be tuned as in Example 14.41 D 
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Example 4.9: Consider the system 

x = sin(x) + 0.1 sin t — 0.015 

with measurement y = x. Letting eo(y) = y, e\{y) = Ay, and e2(y) = 3y, the 
variational equation is 

5x + 3 5x + (4 — cos(x)) Sx + 5x = 

The corresponding observer is illustrated in Figure |H1 d 



STATE ESTIMATES 

40 | , 1 r 




1 1 1 1 1 1 1 

O 5 10 15 20 25 

TIME(SEC) 



Figure 8: x = sin(x) + 0.1 sin t - 0.015 , with y = x. 



4.3 Adding an indifferent system 

The analysis may be further simplified by using superposition to compare the 
system to one whose contraction behavior is known analytically. We illustrate 
this idea on second-order systems using an indifferent added dynamics. In 
principle, the approach can be extended to higher-order systems as well as 
other types of added dynamics. 

Consider the indifferent system ^2] 

5x = ifWx 
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with real and invertible f2(x, x, t). The above corresponds to the second-order 
dynamics 

which is thus itself indifferent. We can write the reduced form of Theorem El 

as 

d f 5± \ f -6(x,x,i) -OQ(x,x,i) \f S± 
Jt\5x ) ~ { I )\5x 

tin- 1 f -srn^-e o\\ ( s± 



I J \ J J \5x 

which thus corresponds to the superposition of an indifferent system with a 
semi-contracting system of rate — © — f2f2 _1 . 

Theorem 4 The reduced form 

<5x + 0(x, x, t)5± + QQ(±, x, t) 5x = 

is semi-contracting with rate —Q — QQ" 1 . The corresponding unreduced form 
(see Theorem^) has an additional contraction rate — ~. 

Example 4.10: Consider the restricted Three-Body Problem [2^ 

1 \ / 1 
-1 x+1 0| x- W(x) 
00/ 



with state x = (x, y, z) T , potential energy V(x) = —v^Jy 2 + z 2 + (—1 + x + y) 2 

(1 — v) \J 'y 2 + z 2 + (x + y) 2 , v the ratio of the smaller mass to the larger mass, 
and the third body of zero mass. 

1 

Using F = I —1 I in Theorem |3J the reduced variational dynamics is 







<5x = -AU<5x 



Using Theorem @J the contraction behavior of the three-body problem is the 
superposition of 

• i\/ AV, which is indifferent for AV > and unstable otherwise. 

• — d ^dt V V AV , where a tightening (relaxing) potential force adds semi- 
contracting (diverging) behavior. D 
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5 Hamiltonian system dynamics 



Consider the general n-dimensional Hamiltonian dynamics with sum 
convention over the free index 

& + l j kh i k x h = -H^ h f h - D{x h 

with external forces fh{x l ,t), damping gain D 3 h (x l ) and the Christoffel term 
llh = \ Rmj (^T + ^ ~ l§£) (section 7.2. in |Ej) of the symmetric 
inertia tensor Hih(x m ) and where we use the convention H m3 = H^j- The 
variational dynamics of the above is 

SS p + ^k^i^^h + 2 4,x k 6x h = - d ^ h , fh) 5x l - ^Sx l x h - DUx h 
ox 1 ox 1 ox 1 

Let us now compute with the above the covariant time-derivative (section 
3.6. in ^Hj and generalized contraction with metric in [12]) of the covariant 
velocity variation 5x J + jj h x 8x l with respect to the metric Hih as 

j t (5x j + ii h x h 5x l ) + ^ h x h (Sx* + l^xHx 1 ) = 

5x J + ^x k x h 6x l + i 3 lh x h 5x l + il h x h 5x l + ^ h x h {5x l + f kl x k 6x l ) = 

Kf kh x k x h 5x l - W h - 7 k hl f k ) 5x l - D{ (5x h + ^xHx 1 ) (10) 

with curavture tensor Kf kh = - + 7^7^ - 7^7^ (section 7.3. in 
[To]), covariant derivative of the external forces ^ — jj^fk, and where we 

have assumed the covariant derivative of D J h , that is + yj k D k — Jih^i 
(see section 3.6 in ^3]), to vanish which is usually the case for mechanical 
damping. 

Note that according section 7.3 in ^3] corresponds K l lkh H jiY k Y h X 3 X 1 to 
the Riemannian or Gaussian curvature of the 2-D subspace span by X and 
Y. Hence the convexity of H orthogonal to x acts as a spring, whose gain 
increases linearly with velocity. 

Also note that (|10[) can be used directly in combination with Theorem 2 
in [TJ] to show under which condition the covariant Hessian of the action (p 
becomes convex, which implies contraction behavior. 
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However we use here Theorem |H] for F = D, which is more general than 
the above (i.e. it allows to assess complex solutions of the Hessian dynamics 
in Theorem 2 in [E]), to conclude: 

Theorem 5 Consider for t > the n- dimensional (n > 1 ) dynamics with 
sum convention over the free index 

x j + l{ h x k x h = -H jh f h - D{x h (11) 

with external forces p{x , t), damping gain D J h (x l ) with covariant derivative 

8 -^+ 1 l k D k h - 1 lDi = 0, and the Christoffel term ^ kh = + ^ - f^) 

of the symetric, u.p.d. and invertible inertia tensor Hih(x m ). 

Let us define a constant X > as the largest singular value Wt > as 

-K ipkh x x + — - % p f h + H ih 4 I ii I -Kji kh x x +-gj- Ijhh + tijh 4 J <X H pl 

(12) 

uritfi curavture tensor K jlkh = H jo ^t ~ + 7&7« ~ 7k7jU- 
VLe can i/ien conclude on contraction rate — ^ + X. 

Also note that taking the double integral of the dynamics leads to an 
exponential Lyaponov energy stability proof for the autonomous case. In 
this sense Theorem 5 represents a variational extension of the classical energy 
based Lyaponov proofs for autonomous systems. Note that at a variational 
energy approach a tightening spring (positive Q) leads to semi-contraction 
behavior. 



Example 5.1 : Let us now illustrate the simplicity of the stability results using 
the inertia tensor as metric. Consider the Euler dynamics of a rigid body, with 
Euler angles x = (ip, 6, (j)) T and measured rotation vector uj in body coordinates 



— sin 6 

x — | i) rose cos^sin^ ) oj (13) 
cos 9 cos ijj 
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The underlying energy is 



h 





After a straightforward but tedious calculation we can compute 

Thus, the Euler dynamics (jl^j) is globally indifferent. Note that this can also be 
seen from the quaternion angular dynamics, whose Jacobian is skew-symmetric |2U| . 
□ 



Example 5.2: The convexity of a the inertia tensor can act similarly to a 
stabilizing potential force in the variational Hamiltonian dynamics Con- 
sider a rotating point mass of mass m on a ball with radius R. The Hamiltonian 
energy is 

mR 2 . T ( 1 



h = x _ . o , x 

2 V sm & J 

with latitude <j> and longitude ^ in x = (cf),i()) T '. The curvature tensor can be 
computed e.g. with MAPLE as 

fc fc _ mfl 2 sin 2 </> / -V, 2 U 
KjlkhX X ~ 2 ^ ^ -0 2 

which scales as the inertia tensor with and is negative orthogonal to the 
velocity and indifferent along the velocity. Hence the convexitiy in Theorem 
[5] acts under motion as a stabilizing spring, that lets two moving neighboring 
trajectories oscillate around each other. D 



Example 5.3: Theorem [5] can be used to define observers or tracking con- 
trollers for time-varying Hamiltonian systems. Consider a two-link robot ma- 
nipulator, with kinetic energy 

1 i . . n, / ai + 2a 2 cos q 2 a 2 cos q 2 + a 3 \ / q% \ 

2 ^ 1 q2 ' V a 2 cos 172 + 03 ^3 / \ ?2 / 
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with a\ = m\l1i + I\ + m2(l\ + / 2 2 ) + h, «2 = "i-2Wc2 and 03 = m2L 2 c2 + I2 and 
—71" < QiiQ2 < 7T- Let us assume that g- 7 is measured and define the observer 

^'+7L(^W = - Rhj fh 

k 3 = ^+d{(q h -q h ) 

with the external forces 

_ / g{mil c i +m e l 1 )cosq 1 + gm e l ce cos(qi + q 2 ) + n \ , , Ai ^ 

//l "V SW« «*(«!+ ft) +75 ) + *» * Q) 

and external torques ri, 75. The above is equivalent to (fTTjl 

+ ll h (m k q h = ~H jh (f h - d hk q k ) - d\q h 



where the covariant derivative of a constant scalar di vanishes. The curvature 
tensor can be computed e.g. with MAPLE as 

. k . h _ (aia3-Q3- a 2) Q 2COsg 2 ( -q\ q x q 2 
aia 3 - a A 2 cos 2 q 2 - % \ q\q2 -qf 

The curvature is for aia^ > a| + a| convex (concave) for — f < §2 < f 
(— I > ^2 or ^2 > § ) and accordingly (de)-stabilizes the dynamics when the 
arm is retracted (extended). Let us now compute the covariant derivative of the 
external forces 

|§ + 1% (fh - d hk q k ) = -k; + 7* (fh - d hk q k 



w ; t V, -,1 _ ^2 _ 2 _ 2 _ (02 cos 92+03)02 singg 1 _ 1 _ 1 _ 
WltU 7ll - 7 12 - 7 2 i - 722 - a^-a^os* q 2 -a^ ' 7l2 ~ 721 ~ 7 2 2 ~ 

0302 sing 2 1 2 _ (oi+2a 2 cos 02)02 sinq 2 

1 — 9 2" ana 71 1 — 1 — 5 3 — • 

0,1(13— 0,2 cos z 02— a| 111 0103—0^08^02—03 

The spring gain k l p stabilizes the system, whereas the supporting force fh — dh k q k 
can (de)-stabilize the system proportional to the magnitude of the supporting 
force. Note that fh — dh k q k can (de)-stabilize a curved system is unavoidable 
since here no constant or parallel (force) vectors exist, whose covariant derivative 
vanishes [Td] . 

Computing A from (|12fl then allows with TheoreomElto compute bounds on the 
velocity q 1 and external forces fh — dh k q k for which global contraction behavior 
can be concluded. 

System responses to a control input t% = (sini, cos5i), initial conditions q l (0) = 
(-f,vr) rad, g*(0) = (3,-3) rad/s, g l '(0) = (-§,tt) rad, <f(0) = (-5,5) rad/s 
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and parameters mi = 1 kg, l\ = 1 m, m e = 2 kg, I\ = 0.12 kgm 2 , l& = 0.5 m, 
I e = 0.25 kgm 2 , l ce = 0.6 m, d 3 h = 5 Nms/rad, fci = 5 Nm/rad are illustrated in 
figure I§1 and 1101 The solid lines represent the real plant, and the dashed lines 
the observer estimate. 

The above observers provide a simple alternative to current design methods [see 
e.g., Berghuis and Nijmeyer, 1993; Marino and Tomei, 1995], and guarantees 
local (for bounded velocities and time-varing inputs) exponential convergence. 

Note that Theoreom[S]can also be used to bound the diverging behavior, caused 
by the concave inertia when the robot arm is pointing inwards, of the double 
inverted pendulum, when no damping or stabilizing potential force is applied. 
□ 

Example 5.4: Biology found a solution to the problem that a time- varying 
supporting torque can destabilize a system. 

Recently, there has been considerable interest in analyzing feedback controllers 
for biological motor control systems as combinations of simpler elements, or 
motion primitives. For instance [I] and ^0] stimulate a small number of areas 
(A, B, C, and D) in a frog's spinal cord and measure the resulting torque angle 
relations. 

Force fields seem to add when different areas are stimulated at the same time 
so that 0] and ^2] propose the following biological control inputs 

n 

= WW) 

1=1 




TIME 



Figure 9: Positions of two-link robot 
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Figure 10: Velocities of two-link robot 

where each single torque ki(£)fn{cp) results from the stimulation of area I in 
the spinal cord. With force measurements the above authors did show that the 
covariant derivative + jf p fhi of fu with respect to the inertia tensor of the 
frog's leg or body is uniformly positive definite. 

Likely candidates for ki(t) are positive upper and lower bounded sigmoids and 
pulses and periodic activation patterns. 

Using Theorem or the discussion in Example 15.31 with sufficient damping then 
allows to compute a maximal cp for which exponential convergence to a single 
motion is guaranteed. 

Note that the achievement of tracking control with a proportional gain k[(t) 
rather than an additional supporting force as in Example 15 . 31 has the advantage 
that the supporting force has no impact on the contraction behavior anymore. 
□ 



6 Concluding remarks 

The research in this paper can be extended in several directions, as the 
development suggests. 

Some of the extensions will likely require the combination of the above 
results with a simplifying metric pre-transformation, as mentioned in section 
3.1. In particular, classical transformation ideas in nonlinear control such 
as feedback linearization and flatness [S| typically use linear time-invariant 
target dynamics, while the framework provided in this paper should allow 
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considerably more flexibility. This, combined with the fact that a metric 
transformation such as 8z = @5x need not be integrable (i.e. does not 
require an explicit z to exist), could potentially lead to useful generalisations 
of these methods. 

Acknowledgement The authors are grateful to Yong Zhao for per- 
forming the simulations and for stimulating discussions, and to Wei Wang 
for thoughtful comments and suggestions. 
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